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Abstract 

We develop an efficient algorithm for evaluating divergent perturbation expansions of field 
theories in the bare coupling constant gB for which we possess a finite number L of expansion 
coefficients plus two more informations: The knowledge of the large-order behavior proportional 
to {—a)''k\k^g'^, with a known growth parameter a, and the knowledge of the approach to 
scaling being of the type c + c'/g^, with constants c,c' and a critical exponent of approach lu. 
The latter information leads to an increase in the speed of convergence and a high accuracy of the 
results. The algorithm is applied to the six- and seven-loop expansions for the critical exponents 
of 0(A^)-symmetric (/)'*-theories, and the result for the critical exponent a is compared with the 
recent satellite experiment. 



\0 1 Introduction 

^O ' The field-theoretic approach to critical phenomena provides us with power series expansions for the 
O^ critical exponents of a wide variety of universality classes. For (/)'*-theories with 0{N) symmetry in 



three dimensions, these expansions have been calculated numerically as power series in the renor- 
malized coupling constant up to seven loops for the critical exponents v and t] |i| and up to six loops 
for the exponent uj governing the approach to scaling. In 4 — e dimensions, exact e-expansions are 
available up to five loops for all critical exponents with O(A^) symmetry [g], cubic symmetry, and 
K*" mixtures of these 0. When inserted into the renormalization group equations, these expansions 
are supposed to determine the critical exponents via their values at an infrared-stable fixed point 
g = g* . The latter step is nontrivial since the expansions are divergent and require resummation, 
for which sophisticated methods have been developed, summarized and applied most recently in Q. 
The resummation methods use the information from the known large-order behavior {—a)^k\k^g^ 
of the expansions and analytic mapping techniques to obtain quite accurate results. 

A completely different resummation procedure was developed recently on the basis of variational 
perturbation theory |^] to the expansions in powers of the hare coupling constant, which goes to 
infinity at the critical point. The resulting strong- coupling theory P] was successfully applied in 



three [^, ^ and 4 — e dimensions |10|, and rendered for the first time an interpolation between 



expansions of 4 — e and 2 -|- e-dimensional theories. This method converges as fast as the previous 
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ones, even though it does not take into account the information on the large-order behavior of the 
expansions. Instead, it uses the fact that the power series for the critical exponents approach their 
constant critical value in the form c-|- c'/^^g, where c, d are constants, and uj is the critical exponent 
of the approach to scaling. The results showed that the latter information is just as efficient in 
increasing the speed of convergence as the information on the large-order behavior. 

We may therefore expect that a resummation method which incorporates both informations 
should lead to results with an even higher accuracy, and it is the purpose of this paper to present 
such a method in the form of a simple algorithm. 

2 The Problem 

The development of our resummation algorithm is based on an improvement of the problem formu- 
lated in [0, 1^ and solved via variational perturbation theory 0]. Mathematically, the problem we 
want to solve is the following: Let 

fL{gB)=J2fkg'h (1) 

k=0 

be the first L terms of a divergent asymptotic expansion 



f{9B)=J2fk9% (2) 



fc=0 

of a function /(^s), which possesses a strong-coupling expansion of the type 

oo 

f{9B)=giY.^kgB'''', (3) 

fc=0 

which is assumed to have some finite convergence radius \gB\ > g's^^ • Suppose that the function is 
analytic in the complex (jr^-plane with a cut along the negative real axis, with a discontinuity known 
from instanton calculations |11, 12 1 to have near the tip of the cut the generic form 



disc/(-5B) = 27ri7(a|9B|)~^''e"'/"l^^l. (4) 

Via a dispersion relation, 

f{gB) = ^z— / dge — r-^ , (5) 

2vr« Jo g's + gB 

or a sufficiently subtracted version of it, this discontinuity corresponds to the large-order behavior 
of the expansion coefficients fk 

fk'='^lkl{-a)''k('[l + 0{l/k)]. (6) 



The constant a is given by the inverse action of the radial symmetric solution to the classical field 
equations. The parameter f3 counts the number of zero modes in the fluctuation determinant around 
this solution. The absolute normalization 7 of the large-order behavior requires the calculation of 



the fluctuation determinant |11]. 



As far as the leading strong-coupling coefficient bo is concerned, this problem has been attacked 



before by Parisi |13| using a resummation method based on Borel transformations in combination 
with analytic mapping techniques. However, when applied to the asymptotic expansions of the 
ground state energy of the anharmonic oscillator, his method converged very slowly, too slow to 
lead to reliable critical exponents, where only five to seven expansion coefficients fk are known. The 
reason is that in Parisi's approach, the corrections to the leading power behavior are failing to match 
the true fractional powers of the strong-coupling expansion (^). 

This deficiency was cured by the strong-coupling theory of one of the authors (HK) in Ref . Q , 
and the subsequent application to critical exponents in Refs. ^ ^ |l^, which showed a surprisingly 
rapid convergence. However, that theory did not take advantage of the knowledge of the large-order 
behavior IM), which can lead to an increase in the speed of convergence and thus of the accuracy of 
theoretical values for the critical exponents. This will be done in the present improved resummation 
method. 

3 Borel Methods 

Basis for this method is the development of a more general Borel-like transformation which will 

automatically guarantee the form of the strong-coupling expansion (|3|) for each approximant fhiOB)- 

Let us first recall briefly the important properties of the ordinary Borel transformations: It is a 

function B{t) associated with f{gB) which is defined by the Taylor series 

00 00 ~ 

B{t) = ^B,t'^Y.^_t'- (7) 

fe=0 k=0 

By dividing the expansion coefficients fk by k\, the factorial growth of /^ is reduced to a power 
growth, thus giving B(t) a finite convergence radius. 

An alternative definition of the Borel transform is given by the contour integral 

where the contour C encloses anticlockwise the negative real axis. Indeed, inserting here (|l|) and 
performing the integral we obtain (j^. 



If f{gB) is an analytic function in the sector 

S^/2+5 = {9B \\gB\<R,\ arg(5B)| < ^/2 + 6} (9) 

of a circle, and satisfies the so-called strong asymptotic condition 

L 

\f{gB)-Y.fkg%\<Ag'^+^a'^+\L + l)\, with a, A > 0, (10) 

fc=0 

then B{t) is analytic in S^ , with a finite radius of convergence t < 1/a. The original function /{gs) 
can be recovered from B{t) by the inverse Borel transformation 

/•oo 

fias) = / e-''B(tgB)dt. (11) 

JO 

Obviously, the inverse transformation can only be performed if B{t) is known on the entire positive 
real t-axis. The Taylor series (0) for B{t), however, converges only inside the circle of radius 1/a. 
Before we can do the integral in (pT|), we must therefore perform a suitable analytic continuation of 
(0) 13 ■ This can be done by reexpanding B{t) in powers of the function nit) defined implicitly by 

1 nit) , , 

This function maps the interval [0, cxd] of the t-axis to the interval [0, 1] of the K-plane. By a proper 
choice of a it is possible to make the unit circle free of singularities. Then we may use the reexpansion 
B{t) in powers of K{t) truncated after n^ , 

L 

i?L(t)^ j;r;fcAi^(t), (13) 

fc=0 

as an approximation to BL(t) on the entire positive real t axis. Inserting this into the inverse 
transformation formula ([ll|) , we obtain an approximation f1{gB) for /{gs), which has the same first 
L expansion coefficients as JLigs) and, in addition, the correct large-order behavior ^. 

How can we incorporate the strong-coupling expansion (^) of /{gs) into the approximation 
ftidB)'^ In the Borel transform B{t), the strong-coupling expansion (y) amounts to a large-t expan- 
sion 

Bit) =t^y ^i"^(^^-^) r(fc^ _ s)b,t-'^^. (14) 

fc=0 

This follows directly by inserting (|3|) into (|^) and integrating each term. Here and in the sequel, 
C denotes a path of integration which encloses anticlockwise the negative real axis in the complex 
plane. 



If the series (|3|) has a finite radius of convergence, the large-t expansion of B{t) is a divergent 
asymptotic one, because of the factor T{kuj — s) in the k-th expansion coefficient. 

It should be stressed, that the relation between the coefficients of the strong-coupling expansion 
(^) and the coefficients of expansion (|I^ is not generally invertible, because of the factor sin Tr{kLO — s) 
which causes the coefficients of negative integer powers of t to vanish. 

Note that, in general, an expansion in the Borel-plane with a power sequence in t as in (^) is 
not sufficient to ensure an expansion in the same powers in the g^-plane as in (^, because of the 
appearance of extra integer powers in gs- This is illustrated by the simple function B{t) = (1 + ty, 
which possesses a strong-coupling expansion in the powers f*~ . If s is non-integer the expansion of 
the corresponding function /{gs) reads 

figs) = J^ dt e-\l + gsty = e'/^-T{s + 1)9% + e^^- ^ —l-±—g-^'^-\ (15) 

and expanding the exponential we see that the sum contains integer powers which are not contained 
in the strong-coupling expansion of B{t). 



It is advantageous to perform a further analytic continuation of the reexpansion (13) which 
enforces automatically the leading power behavior t* of B{t). For this we change (^) to 

L 
i?L(t)^[l-«(t)]-P^J]/lfc/.'=(t). (16) 

fc=0 

The coefficients h^ are determined by using (|l2|) to expand K{t) in powers of t, inserting this into 
(|T6[), reexpanding in powers of t, and comparing the final coefficients with those in (|^). When the 
approximation (|lq) is inserted into (|ll|), we obtain /£(5_b) with the correct leading power behavior 
5(| for large gs- 

Unfortunately, the simple prefactor does not produce the correct subleading powers {gsY' 
of the strong-coupling expansion (^), and we have not been able to find another simple analytic 
continuation of B{t) which would achieve this. 

4 Hyper-Borel Transformation 

A solution of this problem is, however, possible with the help of a generalization of the Borel-Leroy 



transformation to what we shall call a hyper-Borel transformation [15| 



B{y)=Y,Bky\ (17) 

fc=0 

with coefficients 

~ _ r(fc(l/a;-l) + /?o) 
r (k/uj - s/uj) T{(3o) 



4.1 General Properties 

The inverse transformation is given by the double integral 

"y^(l-a;)/a; 



'o y 



;(,,). nM£..,-«.J 



9B 



yt(l-'^)/'^ 



exp 



9B 



Biy), (19) 



as can easily shown with the help of the integral representation of the inverse Gamma function 

-^ = -^ /" dteH-'. (20) 

The transformation possesses a free parameter /?o which will be used to optimize the approximation 
fiigB) at each order L. The power lo of the strong-coupling expansion is assumed to lie in the 
interval < lj < 1, as it does in the upcoming physical applications. 

The hyper-Borel transformation has the desired property of allowing for a resummation of fLigs) 
with the full sequence of powers of gs in the strong-coupling expansion (|3|). To show this we first 
observe that as in the ordinary Borel transform (^, the large-argument behavior of the Gamma 
function known from Stirling's formula 

T{pk + q) '^ =°° V2^'^'^pi-^/^k-'^/'^+'^-P/'^pP''{k\)P[l + 0{l/k)], (21) 

removes the factorial growth (y) from the expansion coefficients fk, and leads to a simple power 
behavior of the coefficients 5^: 



f^ fe— >oo , 

ijfc = const. X 



atj(l-w)^/'^-^l k'^+'^''+'^/'^+'/^[l + 0{l/k)]. (22) 



Thus our transform B{y) shares with the ordinary Borel transform B{t) the property of being 
analytic at the origin. Its radius of convergence is determined by the singularity on the negative real 
axis at 

4.2 Resummation Procedure 

A resummation procedure can now be set up on the basis of the transform B(y) as before. The 
inverse transformation ( |l9|) contains an integral over the entire positive axis, requiring again an 
analytic continuation of the Taylor expansion of B{y) beyond the convergence radius. 

The reason for introducing the transform B{y) was to allow us to reproduce the complete power 
sequence in the strong-coupling expansion (|3|), with a leading power g^ and a subleading sequence 
of powers 5^~ '^, A; = 1, 2, 3, . . . . This is achieved by removing a factor e^f^^ with p,a > from the 



truncated series (|l8|) of our transform B{y). Furthermore by removing a second simple prefactor of 
the form (1 + <jy)^ we weaken the leading singularity in the hyper-Borel complex y-plane, which 
determines the large order behavior (pi). The remaining series has still a finite radius of convergence. 
To achieve convergence on the entire positive y axis for which we must do the integral (|l^), we 
reexpand the remaining series of y in powers of n{y) which is related to y by an equation like 
Eq. (|l^). For simplicity we choose the parameter p = 1, i.e. 

1 K{y) 



y 



(24) 



cr 1 - K{y) 

which maps a shifted right half of the complex y-plane with Re[y] > —l/2a onto the unit circle in 
the complex K-plane. Thus we reexpand B{y) in the following way: 

k=0 



(1 + ^y) 



k+S' 



(25) 



A:=0 fc=0 

The inverse hyper-Borel transform of -B(y) is now found by forming the integrals of the expansion 
functions in (p5|) 

1 / 1 n <^ 



ln{g 



B 



27ri 



dteH 



t+-/3o 



C 



dy 
y 



9B 



exp 



93 



-pcjy _ 



(i + ^y) 



n+S '■ 



(26) 



so that the approximants fligs) may be written as 



11(93) =^hnln{93)- 



(27) 



n=0 



The same functions In{g3) may be used as basis functions for a wide variety of divergent trun- 
cated perturbation expansions fhigs)- The complete list of parameters on which they depend reads 
as follows 

In{93) = In{gB,^,s,p,a,5,Po) = In{(7gB,^, s, p,l, 6, (3o), (28) 

but in the following we shall mostly use the shorter notation Inigs)- The integral representation of 
Inigs) breaks down at s = n, requiring an analytical continuation. For the upcoming applications in 
the large-^B regime it will be sufficient to perform this continuation only in the convergent strong- 
coupling expansion of Inigs)- This is obtained by performing a Taylor series expansion of the 
exponential function in (p6|), which is an expansion in powers of l/^lg. After integrating over t and 



y using (20), we obtain an expansion 



^(5b)=5|E^™'^ 



(29) 



fc=0 



which has indeed the same power sequence as the strong-couphng expansion 
figs) to be resummed. 

The expansion coefficients are 

An) _ (-1)*^ a'-'^^TiPo) 



of the function 



k\ T[{uj - l)k + Po + {l/uJ - l)s] ^ 



k ' 



(30) 



where r^ denotes the integral 



,-(") 



Jo 



S~n kw+n—s—l 



(31) 



This integral is seen to coincide with the Kummer function 

1 /"OO 



U{a,-/,z) 



r(a) Jo 



dye-^V"i(l + y)^-°-\ 



(32) 



in terms of which we can write 



An) 



j, = T{ku! + n — s) U{kLO + n — s,kuj — s — 6 + 1, p). 



(33) 



The latter expression is useful since in some applications the integral (^) may diverge, and re- 
quires an analytic continuation by deforming the contour of integration. Such deformations are 
automatically supplied by choosing other representations for the Kummer function, for instance 



U{a,-f,z) 



vr 



sm7r7 



M(a,7,z) 



r(l + a-7)r(7) 



i_^ M(l + a-7,2-7,z) 
' r(a)r(2-7) 



(34) 



where M{a, 7, z) is the confluent hyper geometric function with a Taylor expansion 

The alternative expression ( 33 ) for ij^ , with (^) and ( p5|) , is useful for resumming various asymp- 
totic expansions, for example that of the ground state energy of the anharmonic oscillator, in which 



case the leading strong-coupling power s has the value 1/3. There the integral representation (31) 
would have to be evaluated for values n = 0, k = 0, where the integral does not exist, whereas 



formula (33) with (34) and (|35|) is well-defined. 

For large k, the integral on the right-hand side of (|3l| ) can be estimated with the help of the 
saddle point approximation. The saddle point lies at 

kuj 
P 



Vs 



(36) 



leading to the asymptotic estimate 



1") '=°° (—) " /" dye-^V'+""'"' [1 + 0(1A)] 



—^ /5-^'^-"+T(A:u; + n-s)[l + 0{l/k)] . (37) 



(n) 

The behavior of the strong-couphng coefficients bj^ for large k is obtained with the help of the 



identity 



T{z)T{l-z) = -^^ (38) 

smvrz 



and Stirling's formula (|2l|), yielding 

k 

^W fc^ioo ^sinvr[A:(a;-l)+/?o + (l/a;-l)s] 



(l-a;)(i-^) 



A;^i[l + 0(1/A:)]. (39) 



The values of the real constants 7, 71 will not be needed in the upcoming discussions, and are 
therefore not calculated explicitly. 

Equation ( |39|) shows that the strong-coupling expansion (|^) has a convergence radius 

(pa)" , .^, 

which means that the basis functions InigB)^ and certainly also /{qb) itself, possess additional 
singularities beside qb = 0. The parameter p will be optimally adjusted to match the positions of 
these singularities. 

4.3 Taylor Series of Basis Functions 

For reexpanding fLigs) in terms of the basis functions Inigs), we must know their Taylor series. 
These are obtained by substituting into ( p6| ) the variable y by qbu' , and expanding the integrand of 
( p6|) in powers of gB- After performing the integrals over y' and t, we find 

00 

In{gB) = Y.ft^9'B, (41) 

k=n 

with the coefficients 

(„) ^ 1 T{Po)T{k/u - sH "^ f-S - n\ {-pf-^-^ , 
^^ u:T{k{l/u-l)+l3^)j^^\ J ){k-n-j)\ ' ^ ' 

The coefficients in the last sum arise from the i-integral: 

^V ] ){k-n-j)\ T{k-n+i)T(n + 5) ]„ 



For large k, the integral may be evaluated with the help of the saddle-point approximation. Using 
this and Stirling's formula (^), we find 

^V J j{k-n-3)\ T{5 + n) 

Inserting this into (H^) and using once more Stirling's formula, we obtain for the expansion coefficients 

In) 

ff, the following factorial growth 



An) k~^oo ( l)"ePr(/3o) ^l/2-f3o^l3o-l+s/ujjJ-l3o+n-3/2-s/uj 



a 



-\k 
k\ 



x[l + 0(lA)]. (45) 

For an optimal reexpansion (^), we shall choose the free parameters of the basis functions 
InigB-,^-, s, p, a, 6, /5o) to match the large-order behavior of the coefficients fk in (^). 

4.4 Convergence Properties of Resummed Series 

We shall now discuss the speed of convergence of the resummation procedure. For this it will be 
sufficient to estimate the convergence of the strong-coupling coefficients 6^ of the approximations 
/lIs'b) against the true strong-coupling coefficients bk in (y). The convergence for arbitrary values of 
qb will always be better than that. Such an estimate is possible by looking at the large-n behavior of 
the expansion coefficients bj^ in the strong-coupling expansion oiIn{gB)- in (p9[). This is determined 
by the saddle point approximation to the integral ij^ in Eq. (|3l]), which we rewrite as 

, . /-oo 

if^= dye-^S'-"l'^(i+i/2^)(l + y)-V^-^-i. (46) 

Jo 

The saddle point lies at 

y, = y^[l + 0(l/V^)]. (47) 

At this point, the total exponent in the integrand is 

-pVs-n In ( 1 + 1 ) = -2^^ [1 + 0(1/ V^)] , (48) 

V VsJ 

implying the large-n behavior 

The strong-coupling coefficients 6^ of the approximations f1{gB) are linear combinations of the 
coefficients b^ of the basis functions InidB)- 

L 



b'k=E€'hn. (50) 



n=0 

10 



The speed of convergence with which the 6^'s approach b^ as the number L goes to infinity is governed 
by the growth with n of the reexpansion coefficients /i„ and of the coefficients bj^ in Eq. (|49| ) . We 
shall see that for the series to be resummed, the reexpansion coefficients /i„ will grow at most like 
some power rf , implying that the approximations 6^ approach their L ^ oo -limit b^ with an error 
proportional to 

The leading exponential falloff of the error e~'^^f^ is independent of the other parameters in the 
basis functions /„(gB,u;,p, p, o", (^, /3o) which still need adjustment. This is the important advantage 
of the present resummation method with respect to variational perturbation theory P, H] where the 
error decreases merely like e"™"^*^*^ '^ with 1 — w close to 1/4. 

The nonexponential prefactor in Eq. (^) depends on the parameters in In{gB,^,P, P^f^^^^Pa)- 
Some of them are related to observables, others are free and may be chosen to optimize the conver- 
gence. 

4.4.1 Parameters s and uj 

The perturbation expansions for the critical exponents are power series in the bare coupling con- 
stant qb whose strong-coupling limit is a constant |8|, ^. The same is true for the series expressing 
the renormalized coupling constant g in powers of the bare coupling constant. This implies that 
the growth parameter s for the basis functions Inigs) is equal to zero in all cases. The constant 
asymptotic values are approached with the subleading powers l/g^~^ , where w is a universal ex- 
perimentally measurable critical exponent. 

4.4.2 Parameter a 

In the ordinary Borel-transformation, the parameter a in the large-order behavior of the expansion 
coefficients fk in Eq. (0) , which is determined directly by the inverse value of the reduced action of the 
classical solution to the field equations, specifies also the position of the singularity on the negative 
t axis in B{t). In our transform B{y), the singularity position of the singularity is proportional to 
a, with an a;-dependent prefactor. It lies at [see Eq. (|23|)] 

a = auj{l - uoy/'^-K (52) 



This value of a ensures that the expansion coefficients /^ of the basis functions Inigs) in Eq. (^ 
grow for large k with the same factor (—a) as the expansion coefficients for /{gs) in Eq. (P). 



11 



The conformal mapping ( p4| ) maps the singularity at t = — l/cr to k = oo, and converts the cut 
along the negative into a cut in the K-plane from 1 to cxd. The growth of the reexpansion coefficients 
hn with n is therefore determined by the nature of the singularity of B{y) at oo. 

In the upcoming applications to critical exponents it will turn out that the value ( p2| ) following 
from the inverse action of the solution to the classical field equations and uo will not yield the fastest 
convergence of the approximations f^igs) towards f{gB), but that a slightly smaller value gives 
better results. This seems to be due to the fact that the classical solution gives only the nearest 
singularity in the hyper-Borel transform B{y) of /{gs)- In reality, there are many additional cuts 
from other fluctuating field configurations which determine the size of the expansion coefficients 
fk at pre-asymptotic orders k. Since the few known /^'s are always pre-asymptotic, they are best 
accounted for by an effective shift of the position of the singularity into the direction of the additional 
cuts at larger negative y, corresponding to a smaller a. 

4.4.3 Parameter p 



According to Eq. (|40|) , the parameter p determines the radius of convergence of the strong-coupling 
expansion of the basis functions In{gB)- It should therefore be adjusted to fit optimally the corre- 
sponding radius of the original function /{gs)- Since we do not know this radius, this adjustment 
will be done phenomenologically by varying p to optimize the speed of convergence. Specifically, 
we shall search at each order L for a vanishing highest reexpansion coefficient hi or, if it does not 
vanish anywhere, for a vanishing derivative with respect to p: 

H{P) = 0, or ^ = 0. (53) 

dp 

AAA Parameter 5 

From Eq. (45)] we see that the parameter 5 influences the power k^ in the large-order behavior (|6|). 
By comparing the two equations, we identify the growth parameter (3 of Inigs) as being 

f3 = S- Pq- 3/2- s/io + n. (54) 

At first it appears to be impossible to give all basis functions Inigs) the same growth power P 
in (|45|) by simply letting 6 depend on the order n as required by (H). If we were to do this, we 
would have to assign to 5 the value 

5 = 6n = P + (3o + 3/2 + s/uj- n, (55) 



12 



which depends on the index n of the function InidB)-, and this means that we perform an analytical 
continuation of the power series expansion of B{y) by reexpanding it as follows: 

oo oo 

B{y) = Y^ Bkv' = e-P^y{l + ay)-' J^ /^fc(ay)^ (56) 

fe=0 fc=0 

But the series in this formula which is obtained from the series of B{y) by removing a simple factor 
still has the same finite radius of convergence and could not be used to estimate B{y) for large values 
of y needed to perform the back-transform ([l9|). It is, however, possible to sidetrack this problem 
by letting p grow linearly with the order L. Then the exponential factor of ( |56D suppresses the 
integrals over y for large y sufficiently to make the divergence of the reexpanded series (56) at large 



y irrelevant. If we determine p from the condition (^), the growth of p with L turns out to emerge 
by itself. 

4.4.5 Parameter Pq 

The parameter /Sq has two effects. From Eq. ( ^0| ) we see that for 

, > ,^ , '^ + W"-')' (57) 

1 — UJ 
the signs of the strong-coupling expansion coefficients start to alternate irregularly. This irregularity 
weakens the convergence of the higher strong-coupling coefficients 5^ with k > kc against bk- The 
convergence can therefore be improved by choosing a /?o which grows proportionally to the order L 
of the approximation. 



In addition, Pq appears in the power of A; in (45), which is a consequence of the fact that it 
determines the nature of the cut in B{y) in the complex y-plane starting at y = —1/a [see Eq. (25)]. 

If we expand both sides of ( pSJ ) in powers of k = cry /{I + ay) and compare the coefficients of 
powers of k, it is easy to write down an explicit formula for the reexpansion coefficients hn in terms 
of the coefficients Bj of B{y) by 

where Bj are obtained from the original expansion coefficients fk of /{gs) by relation ([l8|). 

Before beginning with the resummation of the perturbation expansions for the critical exponents 
of 0^-field theories, it will be useful to obtain a feeling for the quality of the above-developed resum- 
mation procedure, in particular for the significance of the parameters upon the speed of convergence. 
We do this by resumming the often- used example of an asymptotic series, the perturbation expansion 
of the ground state energy of the anharmonic oscillator. 

13 



4.5 Resummation of Ground State Energy of Anharmonic Oscillator 

Consider the one-dimensional anharmonic oscihator with the Hamiltonian 

H = ^ + m'^+gBx\ (59) 

In this quantum mechanical system, there is no need to distinguish bare and renormalized coupling 
constants, but since the previous resummation formulas were all formulated in terms of qb we shall 
keep this notation also here. The ground state energy has a perturbation expansion 

oo 

E^''H9B) = Y.fk9B, (60) 



whose coefficients can be calculated via the Bender- Wu recursion relation |16] to arbitrarily high 
orders, with a large-order behavior 

A = -/5 A;!(-3)^fc-V2[i + o(i/fc)]. (61) 

By comparison with ^ we identify the growth parameters 

a = 3, (5 = -1/2. (62) 

A scale transformation x — > g^'^x applied to the Hamiltonian ( p9| ) reveals the scaling property |jl^ 
for the energy as a function of gs and m^: 

E{m',gB) = g'J'Eig-'^^m^, 1). (63) 



Combining this with the knowledge |17] that E{m , 1) is an analytic function at m = 0, we see that 



E{l,gB) possesses a power series expansion of the form (^), with the parameters 

s = 1/3, uj = 2/3. (64) 



Inserting the latter number together with a from Eq. (62) into (52), we identify 



The ground state energy E^^>{gB) obeys a once-subtracted dispersion relation |17|: 

2 -^ Jo 9b 9b + 9B 



The perturbation expansion (|60|) is obtained from this by expanding 1/(5^ + ^b) in powers of gB, 
and performing the integral term by term. This shows explicitly that the large-order behavior ( |6l|) 
is caused by an imaginary part 



ImS(0)(_|^^|)^ /e; / 1 ^„i/3|g^l ^ 
V 71" V 3\gB\ 
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near the tip of the left-hand cut in the complex (j^-plane, in agreement with the general form (|5| 
associated with the large-order behavior (P). 

Let us now specify the parameter 5. We shall do this here in an n-dependent way using Eq. (|55[) 



which now reads with (64): 

S = 6n = Po + 3/2 - n. (68) 

The corresponding basis functions 

InigB, 2/3, 1/3, p, 2/V3, /3o + 3/2 - n, /3o), (69) 

have then all the same large-order growth parameter /3 in (|^). 

The two parameters p and Pq are still arbitrary. The first is determined by an order-dependent 
optimization of the approximations via the conditions (^). The best choice of /3o will be made 
differently depending on the regions of gs- 

Let us test the convergence of our algorithm at small negative coupling constants qb, i-c, near 
the tip of the left-hand cut in the complex ^r^-plane. We do this by calculating the prefactor 7 
in the large-order behavior (|6|). In this case the convergence turns out to be fastest by giving the 
parameter /?o a small value, i.e. /?o = 2. With the large-order behavior (EBh of the basis functions 



Inigs), we find the resummed functions fLigs) of Lth order 'Yln=Q^nIn{gB) to have a large-order 

behavior ^ with a prefactor 

e'T(/?o 



IL 



v/2vrr(5) ^^Q 

The values of these sums for increasing L are shown in Fig. |l|. They converge exponentially fast 
against the exact limiting value 

7 = \/^' (71) 

with superimposed oscillations. The oscillations are of the same kind as those observed in variational 
perturbation theory for the convergence of the approximations to the strong-coupling coefficients hk 
(see Figs. 5.19 and 5.20 in Ref. |^]) Also here, the strong-coupling coefficients 6^ converge exponen- 
tially fast towards bk, but with a larger power of L in the exponent of the last term ~ g-constxVL jggg 
Eq. (51)], rather than « g-constxL £qj, variational perturbation theory [see Eq. (5.199) in Ref. Q]. 



This is seen on the right-hand side of Fig. [l|. 

We have applied our resummation method to the first 10 strong-coupling coefficients using the 
expansion coefficients fk up to order 70. The results are shown in Table 1. Comparison with a 
similar table in Refs. [|18| , |5[ shows that the new resummation method yields in 70th order the same 
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n 


bn 





0.667 986 259 155 777 108 270 962 02 


1 


0.143 668 783 380 864 910 020 319 


2 


-0.008 627 565 680 802 279 127 963 


3 


0.000 818 208 905 756 349 542 41 


4 


-0.000 082 429 217 130 077 219 91 


5 


0.000 008 069 494 235 040 964 75 


6 


-0.000 000 727 977 005 945 772 63 


7 


0.000 000 056 145 997 222 351 17 


8 


-0.000 000 002 949 562 732 709 36 


9 


-0.000 000 000 064 215 331 956 97 


10 


0.000 000 000 048 214 263 789 07 



Table 1: Strong-coupling coefficients h^ of the 70-th order approximants EjQ{g) = Yln=o^nIn{g) to 
the ground state energy E^{g) of the anharmonic oscillator. They have the same accuracy as the 
variational perturbation-theoretic calculations up to order 251 in Refs. |15, ^]. 

accuracy as variational perturbation theory did in 251st order. In all cases the optimal parameter p 
turns out to be a slowly growing function with L. 

In the strong-coupling regime, the convergence is fastest by choosing for (3q an L-dependent value 



Po = L. 



(72) 



Note that this choice of (3q ruins the convergence to the imaginary part for small negative qb which 
was resummed best with /3o = 2. 

4.6 Resummation of Critical Exponents 

Having convinced ourselves of the fast convergence of our new resummation method, let us now turn 
to the perturbation expansions of the 0(A^)-symmetric (/)^-theories in powers of the bare coupling 
constant qb, defined by the euclidean action 



A = Jd^^xS^^iOMx)? + lmlcj>l{x) + 2n9B[4{x)f 



(73) 



in D = 3 dimensions. The field (j)o is an A^-component vector (j)^ = {(J)q,(J)q, . . . , (J)q), and the action is 
0(A^)-symmetric. We define renormalized mass m and field strength by parametrizing the behavior 
of the connected two point function G^'^' in momentum space near zero momentum as 



G^^\p,a;-p,f3) = Z^- 



^al3 



(74) 



The renormalized coupling constant g is defined by the value of the connected four-point function 
at zero momenta: 



G(^)(0,a;0,/3;0,7;0,5) = m "^ ^Zlg{6ap6^s + ^a-ySps + SaS^p^)- 



(75) 



16 



If we introduce the dimensionless bare coupling constant g^ = gs/m, the critical exponents are 
defined by 



v{9b) =9b-, — fog^,/., 
dgB 



T^{9- 



B) 



9B 



;^— log— ^. 
dgs rn^ 



ml 
,2- 



(76) 



The following expansions for the critical indices in the bare dimensionless coupling constant are 
available [^] in the literature for all O(A^): 

rj^gs) = (16/27 + 8N/27)gl + (-9.086537459 - 5.679085912iV - 0.5679085912A^2)^| 
+(127.4916153 + 94.77320534iV + 17.13477557V2 + 0.8105383221 A^^)^!^ 
+(-1843.49199 - 1576.46676iV - 395.2678358A^^ - 36.00660242iV3 

-1.026437849A^'')5| 
+(28108.60398 + 26995. 87962Af + 8461.481806A^^ + 1116.246863iV^ 

+62.8879068iV^ + 1.218861532A^^)5| , (77) 

2 - i^'H9b) = <7b(2 + iV) + (523/27 + 316iV/27 + N^)gl 

+(229.3744544 + 162.8474234iV + 26.08009809Af^ + A^^)g| 

+(-3090996037 - 2520.8487517V - 572.3282893iV^ - 44.32646141Af^ - N'^)g% 

+(45970.71839 + 42170.32707Af + 12152. 70675 A^^ + 1408.064008iV3 

+65.97630108Af'^ + N^)9% 
+(-740843.1985 - 751333.064A^ - 258945.0037iV2 - 39575.57037A^^ 

-2842.8966A^^ - 90.7145582iV^ - N^)g% . (78) 

In addition, seventh order coefficients have been calculated for A^ = 0, 1, 2, 3: Q: 



.(7) 



' -45387.48927 " 








' -12792269.773 " 




' Af = ■ 


-114574.4876 
-241424.7646 


> 9b, 


V 


-1 (7) = , 


-33711416.972 
-73780809.849 


> g\ for < 


N = \ 

N = 2 


-454761.4731 








-143831857.01 




N = 3 



(79) 



When approaching the critical point, the renormalized mass m tends to zero, so that the problem 
is to find the strong-coupling limit of these expansions. In order to have the critical exponents 
approach a constant value, the power s in Eq. (0) must be set equal to zero. 

In contrast to the quantum-mechanical discussion in the last section, the exponent a; governing 
the approach to the scaling limit is now unknown, and must also be determined from the available 
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perturbation expansions. As in Ref. [^, §, we solve this problem by using the fact that the existence 
of a critical point implies the renormalized coupling constant g in powers of qb to converge against 
a constant renormalized coupling g* for m ^ 0. The expansion of g{gB) is known up to six loops 
[||] for all 0(iV): 

gigs) =gB + (-8 - N)gl + (2IO8/27 + 514/27iV + N^)gl 

+(-878.7937193 - 312.63444671iV - 32.54841303iV2 - N^)g% 
+(11068.06183 + 5100.403285iV + 786.3665699iV2 + 48.21386744iV^ + iV'^)^! 
+(-153102.85023 - 85611. 91996iV - 17317.7025iV2 - 1585.114189iV^ 

-65.820362037V^ - N^)gB 
+(2297647.148 + 1495703.313iV + 371103. 0896iV2 + 44914.04818iV3 

+2797.291579A^^ + 85.21310501iV^ + N^)g'^B ■ (80) 

The convergence against a fixed coupling g* occurs only for the correct value of uo in the resummation 
functions In{gBi^, s, p, a, 6, /?o). At different values, g{gB) has some strong-coupling power behavior 
g^g with s 7^ 0. We may therefore determine a; by forming from ( |80|) a series for the power s, 



s = —T\ = — g [gB), 81 

dloggB g 

resumming this for various values of to in the basis functions, and finding the critical exponent uj 

from the zero of s. Alternatively, since g{gB) ^ 5*) we can just as well resum the series for —gs, 

which coincides with the /3-function of renormalization group theory [not to be confused with the 

growth parameter /3 in (^)] 

^(gs) = -gs'-^. (82) 

dgB 

If we denote its strong-coupling limit by (3* , 

P* ^ (3{gB)\gs^oo, (83) 

we resum the expansion for l3{gB) to form the approximations 

L 

higs) = Y^ KInigB,uj), (84) 

n=0 

and plot the strong-coupling limits of the Lth approximations /?£ for various values of uj. This is 
shown in Fig. g. From these plots we extract the critical exponent m by finding the w-value for which 
the approximations /J^ extrapolate best to zero for L — > oo, taking into account that the convergence 
is exponentially fast with superimposed oscillations. These w-values are called lol- 
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For these resummations, we must of course specify the remaining parameters in the basis func- 
tions In{gB,uj, 0, p, a, 6, (3q). This can, in principle, proceed as in the case of the anharmonic oscillator. 
The parameter a is determined from the action of a classical instanton solution (pdx) of the field 



equations, and has for all expansions the A^-independent value |11] 



327r 
a = -— = 1.32997 , (85) 

h 

where Ip = f d^ x[(/)c{x)]'p are integrals over powers of (pd^). 

To determine the parameter 5, we recall the remaining growth parameters /? and 7 of the large- 
order behavior (|6|) of the perturbative series for the critical exponents. The growth parameter /3 is 
given by the number of zero modes in the fluctuation spectrum around this classical solution: 

Pr, ] ( 3 + N/2 ] 

/?,-! \ = \ A + N/2 \ (86) 

Pp J [ 4 + iV/2 J 

The prefactors 7 in (|6|) requires the calculation of the fluctuation determinants around the classical 
solution, which yields in the case of the /^-function 

^^- r(iV/2 + 2) UJ U" J ' ^ ■ ^ ^ 

where D^ and Dt are characteristic quantities of the longitudinal and transverse parts of the fluc- 
tuation determinant, respectively. Their numerical values are ||11| 



Dl 


Dt 


h 


h 


h 


^3 


10.544 ± 0.004 


1.4571 ± 0.0001 


31.691522 


75.589005 


659.868352 


13.563312 



The parameters 7^, 7,^-1 are obtained from 7^ by: 






7. = 7/3T7^77^^' 7.-=7/3i^^^(^-l)4vr-|, (89) 



where I2 = (1 — D/A)l4 and H^ is listed in (88). Note that the expansions in powers of the 



renormalized coupling constant g have the same parameters a and /3, but different parameters 7/?. 
These differ from the above a's by a common factor: 

7i? = 7e-(^+«)/". (90) 

From Eq. (85), the parameter a is found using relation (p2|). It turns out, however, that this value 
does not lead to an optimal convergence. This can be understood qualitatively by observing that 
the large-order behavior of the expansion coefficients of the critical exponents and of the /3-function 
in powers of the bare coupling constant qb is not nearly as precocious in reaching the large-order 
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form (^) as the corresponding expansions in powers of the renormahzed couphng constant g (see 
Fig. 1 in Ref. Q). The lack of precocity here is illustrated for the expansion coefficients /3fc of the 
/3- function in Table 2, which gives the ratios of (3^ and their leading asymptotic estimates [5'^: 



N 





1 


2 


3 


k 


Pk/Pr 


Pk/Pr 


Pk/P^' 


Pk/Pr 


2 


0.57 


0.45 


0.35 


0.27 


3 


0.61 


0.45 


0.32 


0.22 


4 


0.73 


0.51 


0.34 


0.22 


5 


0.89 


0.61 


0.40 


0.25 


6 


1.07 


0.73 


0.47 


0.29 


7 


1.26 


0.88 


0.56 


0.34 


7/3 


110.0 


97.0 


75.5 


53.2 



Table 2: First six perturbative coefficients in the expansions of the /3-function in powers of the 
bare coupling constant qb, divided by their asymptotic large-order estimates {—a)^klk^i^. The 
ratios increase quite slowly towards the theoretically predicted normalization constant 7/3 in the 
asymptotic regime. 



Pk/PT = Pk/kK-afkf'. 



(91) 



The first six approach their large-order limits quite slowly. For this reason we prefer to adapt a 
not from a by Eq. (52), but by an optimization of the convergence. Since the reexpanded series 
converges for fixed values of 6 and a it is reasonable to determine these parameters by searching for 
a point of least dependence in largest available order L. This is done by imposing the conditions 



dKL „ , (Pkl „ 
and , o = 



(92) 



da da"^ 

to determine both parameters 5, a, where kl denotes the Lth approximation to any exponent 7, v 
or r]. In accordance with the discussions in section 2.1.2 this procedure provides a value of a which 
is smaller than that given by (|5^). 

After trying out a few choices, we have given the parameters P and p the fixed values 1 and 10, 
respectively, to accelerate the convergence. 

The results for the critical exponents of all 0(A^)-symmetries are shown in Figs. |2|-g and Table 3. 



The total error is indicated in the square brackets. It is deduced from the error of resummation 
of the critical exponent at a fixed value of uj indicated in the parentheses, and from the error Aco of 
a;, using the derivative of the exponent with respect to lo given in curly brackets. Symbolically, the 
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relation between these errors is 



[...] = {...) + Au;{...}. 



(93) 



The accuracy of our results can be judged by comparison with the most accurately measured critical 
exponent a parametrizing the divergence of the specific heat of superfiuid helium at the A-transition 
by \Tc — T|~". By going into a vicinity of the critical temperature with AT ~ 10^* K, a recent 
satellite experiment has provided us with the value |20] 



a 



-0.01285 ± 0.00038. 



(94) 



Our value for v in Table 3 is 



I, = -0.6704 ± 0.007 



(95) 



and yields via the scaling relation a = 2 — 3iy: 



a 



-0.0112 ±0.0021, 



(96) 



in good agreement with the experimental number (|9J). A comparison with other experiments and 
theories is shown in Fig. |^, showing that our result is among the more accurate ones. 

A remark is necessary concerning the errors quoted in this paper. We do not know how to esti- 
mating precisely the errors which can appear in an involved numerical approximation scheme such 
as the one presented here. Our estimates are based on the range of critical exponents which can be 
reached by reasonably modifying the parameters in the calculations. What may be considered as 
reasonable is a somewhat subjective procedure. As such, our error estimates follow the rule of maxi- 
mal optimism, and are probably underestimated. This is, however, not uncommon in resummations 
of divergent power series of critical exponents. 



n 


7 


V 


u 


UJ 



1 

2 

3 


1.1604[8] (4) {0.075} 
1.2403[8] (4) {0.110} 
1.3164[8] (5) {0.033} 
1.3882[10] (7) {0.210} 


0.0285[6] (4) {0.037} 
0.0335[6] (3) {0.043} 
0.0349[8] (5) {0.042} 
0.0350[8] (5) {0.043} 


0.5881[8] (4) {0.075} 
0.6303[8] (4) {0.065} 
0.6704[7] (4) {0.098} 
0.7062[7] (4) {0.110} 


0.803[3] {1} 
0.792[3] {1} 
0.784[3] {1} 
0.783[3] {1} 



Table 3: Critical exponents of the 0(A^)-symmetric cj)'^ -theory from our new resummation method. 
The numbers in square brackets indicate the total errors. They arise form the error of the resum- 
mation at fixed values of uj indicated in parentheses, and the errors coming from the inaccurate 
knowledge of to. The former are estimated from the scattering of the approximants around the 
graphically determined large-L limit, the latter follow from the errors in to and the derivatives of the 
critical exponents with respect to changes of iv indicated in the curly brackets. 
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